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Abstract 



In this paper, the approach to asymptotic electromagnetic fields introduced 
by Goldberg and Kerr (J. Math. Phys. 5 (1964)) is used to study various 
aspects of Lorentz Covariant Gravity. Retarded multipole moments of the 
source, the central objects of this study, are defined, and a sequence of con- 
servation equations for these are derived from the conservation of the energy- 
stress-momentum tensor of the source. The solution of the linearized Einstein 
field equation is obtained in terms of the retarded moments for a general 
bound source field, correct to 0{r~^). This is used to obtain the peeling~off 
of the linearized field, and to study the geometric optics approximation for 
the field and for the energy-momentum of the field, given by the Landau- 
Lifschitz pseudotensor. It is shown that the energy-momentum 4- vector splits 
into the 'total radiated 4-momentum' and the 'bound 4-momentum of the 
source', similar to the case of the electromagnetic field. The role played by 
the conservation equations in studying the radiative behaviour of the field 
is stressed throughout. In addition, in the case of a source which has only 
retarded pole, dipole and quadrupole moments, it is shown how to derive the 
arbitrary dependence of the field on a null coordinate. This allows comparison 
with the solutions for linearized gravity obtained by other authors. 



I. INTRODUCTION 

In previous articles @JTI|, the Goldberg-Kerr (GK) |^ approach to the electromagnetic 
field was used to study various aspects of asymptotic electromagnetic fields due to bounded 
sources. In particular, the geometric optics approximation and the arbitrary dependence of 
the field on a null coordinate were established, and the relation of the multipole moments 
of the field to certain retarded moments of the source (integrals over null hyperplanes) was 
determined. One would like to be able to determine the same relationship in the case of 
the gravitational field. However, the non-linear nature of Einstein's field equations prohibits 



*e-mail: nolanb@ccmail.dcu.ie 



1 



this, and one must resort to approximation techniques to obtain results in this vein. Such 
techniques form the keystone in the study of gravitational waves flp^Of. 

The simplest approximation is linearized General Relativity (GR), used in the analysis 
of weak gravitional waves [§,0. Here, the gravitational field is the Riemann tensor obtain 
from a first order perturbation to Minkowski space-time, quantified by writing the metric 
as 



9ij = Vij+lij- (1-1) 

Tensor indices are raised and lowered with the background Minkowski metric rjij, and jij 
and its derivatives are taken to be small of first order, so that any products among these 
terms which arise are neglected. From the metric ( |1.1D the linearized Riemann and Einstein 
tensors can be calculated, and the ensuing linear field equations tackled, with or without 
a non-zero right hand side. Such discussions are usually augmented by a pseudotensor or 
effective stress tensor description of the energy-momentum of the gravitational field. 

The aim of this paper is to study aspects of linearized GR using the GK approach to 
retarded fields due to bound sources. This involves obtaining the solution to the linearized 
field equations explicitly in terms of certian integrals over the source. In full GR, the linear 
approximation certainly does not hold all the way down to the source, but is usually applied 
in the distant wave zone. Hence, to distinguish these cases, the theory studied here will be 
referred to as Lorentz Covariant Gravity (LCG). However, since most of the results presented 
here are of an asymptotic nature, they will apply to the distant wave zone of relativistic 
sources. 

The fundamental object of this study is •jij, and the field equations that we consider may 
be written 

□ 7\ = -167rT,,, (1.2) 
7"^, = 0, (1.3) 

where 7*^^- is the trace reversed part of 7ij, and the trace reversal of any symmetric second 
order tensor is 

f:, = h-\n^,f, f = v''fv- (1-4) 

□ is the d'Alembertian of fiat space-time, and Tij is the energy-stress-momentum tensor of 
the source, which is conserved; 

T{j = 0. (1.5) 
The gravitational field in this theory is given by 

Lijkl = -{lil,jk + 7jfc,?7 — 7ifc,j7 — 'ljl,ik) , (1-6) 

where jij is obtained from the trace reversal of 7*^^-. 

The field due to an isolated, extended, bound source will be considered. Thus the energy- 
stress-momentum tensor Tij of the source obeys the following conditions 0. 
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(1) There exists a time-like world line C : = z^{t) with unit tangent = dz^/dr, 
where r is proper time along C. 

(2) There exists a scalar function /i(r) such that T*-'(x) = for all points x such that 

{x' - z\T))ix, - z,iT)) > max {o, h^r) - [v^ix' - z\r))f} . (1.7) 

Thus in an instantaneous rest frame (IRF) of C, T*-^ vanishes in the infinite region 
extending to spatial infinity which is bounded by the future and past null cones with apex 
z^{t) G C and the cylinder |x — z(r)| = /i(t). The conditions (1) and (2) above define the 
world tube W which is the support of T*-', and in addition the weak assumption of 'strong 
vanishing on the boundary of is made, that is it is assumed that T^^ and a sufficient 
number of its derivatives vanish on the boundary of W. 

Using the world-line C, coordinates and an associated null tetrad for Minkowski space- 
time Ai can be set up in the following way (see |[T8|fl). Let a;" be any point of Ai. The 
intersection of the past null cone with vertex x"' and the world line C is a unique point, 
and defines a unique value u of r, in terms of which is referred to as retarded time. 
Tensor functions of u defined on C then become tensor fields on Ai by parallel transport 
up the null cone N to e.g. f°(a;) = v°'{T)\r=u- a"' ,b°' , respectively are used to denote 
the second (acceleration), third and fourth retarded time derivatives (denoted by a dot) of 
z°'{u). k"" (null) is defined by 

x" = z\u) + rk" , (1.8) 

where the normalization kav"" = — 1 is chosen. Defining the stereographic coordinates (see 
appendix A of [jTH[ for the definition of the polar angles in this case) 

C = e^<^cot^, (1.9) 

one can write 

r = + CC, 1 - CC, C + C, -^(C - 0) , (1.10) 

where 

p = (1 + ay - (1 - cc)^^' - cy - iv') - cy + iv^) . (i.n) 

One can also give 

9,^9 d „ d ,9 d dk" d 

— = k''^—, — = v'^- hra-fc— , Tr=r^^T—, (1-12) 

or ox"- ou ot" or oQ OC, ox"- 



- - (1 + + ^ 1^7^ + 1^7^ h (1-14) 



dx^ 'du ' ' ' '''dr 2r \ dC dC dC dC 
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and 

^ dx"-dx^ r"^ ~^ i^dr"^ ~^ r dr^ r du dudr ' ^ ^ 

where A = p'^d'^/d(d( is the Laplacian on the unit sphere. 

The retarded distance is r = — t>a(x° — z^-^u)), and the following derivatives obtain: 

u,a = -ka, (1.16a) 
Vb,a = -abka, (1.16b) 
r,a = -Va + {l + ra-k)ka, (1.16c) 

h,a = -{Vab + Vah + ^a^^b " (1 + ra-k)kah) . (1.16d) 
r 

Here, a-k = aak"'. From these, it is seen that k"" is tangent to geodesies, with r an affine 
parameter along its integral curves. One can use (|1.16a|) to differentiate tensor functions of 

u; 

Tt{u),a = -ftka. (1.17) 

The field equation ([L^) is solved by (see Synge [|16|, p. 407) 



7\(x)=4/ T,,da = A f ^fy'^y^ d^y, (1.18) 
JNnw J |x — y| 

where da is the 'absolute two-content of the three cell on the null cone' (Synge, [1^ p. 429) 
which reduces to the familiar form in an IRF of C. (Equations capped with an asterisk hold 
in the IRF, and boldface letters represent spatial co-ordinates x", a = 1, 2, 3.) Events in W 
are labelled ?/* = (?/°, y), and the integration is over the intersection of the past null cone 
with vertex with the world tube W. 

In §2, the basic variables of the study, the retarded multipole moments of the source are 
defined. A sequence of conservation equations for these moments are derived using ( p..5[ ), 
and their space-time derivatives are calculated. The solution 7*^^ is obtained to 0{r~^) for 
a completely general source. 

In 53, the linearized curvature tensor is calculated. This allows verification of Sachs' 



peeling-off theorem for linearized gravity [|1^, but for the more satisfactory case of a field 



due to an extended bounded source, rather than one with its source confined to a time-like 
world line ^ . The role played by the conservation equations of §2 in determining the Petrov 
type of the coefficients of the field is stressed, particularly for the leading order (radiative) 
part of the field. The results of this section are given in part in the Newman-Penrose (NP) 
formalism. 

The geometric optics approximation to different aspects of the field is considered in §4. 
First, it is demonstrated how the leading order term in the curvature tensor may be consid- 
ered an approximate solution to the vacuum field equations which obeys the usual geometric 
conditions of a radiation field (algebraic type N with shear-free twist-free geodesic rays). Sec- 
ond, the Landau-Lifschitz pseudotensor is used to discuss the energy-momentum of the 
field. As in the case of the electromagnetic energy tensor @Jl^, it is shown that the energy- 
momentum pseudotensor may be split into two parts, each separately conserved, rad^ij and 
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remtij 



tij 



^Uj. The first term rad^ij has the form associated with a geometric optics 
field. From it, the 'radiative 4- momentum' rad-P* is constructed. This is associated with a 
'cloud of gravitons' which get radiated away to J'~^, while the corresponding rem-P* is shown 
to determine a bound 4-momentum for the field ITHl. In the course of this discussion the 



Sommerfeld outgoing radiation conditions are discussed, and an interesting interplay be- 
tween these and the first of the series of conservation equations for the retarded moments 
of the source is pointed out. 

As mentioned above, it is possible to relate the multipole moments of the field to the 
retarded moments of the source in the case of electromagnetism. One expects the same 
to be true in the case of LCG, and in §5 a decomposition of the field, in the case of a 
pole-dipole-quadrupole source, is obtained. This decomposition would allow one to obtain 
the stated relation, and also brings out the expected arbitrary dependence of the field on 
a null coordinate 
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The resulting alternative form for the solution 



of CO) allows 



comparison with the solutions obtained by other authors 



II. MATHEMATICAL PRELIMINARIES 

In this section, the retarded multipole moments of the source Tij are defined, and the con- 
servation equations which can be deduced from (|1.5| ) are given. The space-time derivatives 
of the first few moments, which will be of use in the following sections, are derived. 

The Goldberg-Kerr formalism is based on the intuitive notion that as the field point 



moves infinitely far away from the source, the domain of integration in ( |1.18| ) is asymptoti- 
cally planar; the portion of the past null cone with apex x"' which intersects W flattens out. 
Thus the central feature of the GK analysis is the use of retarded multipole moments of the 
current 4-vector, defined by 

%j = J^T,,dJ:, (2.1a) 
%j.,k^-k„ = J^Tij^kr--^k„dJ:, ra = 1,2,3,... (2.1b) 



Eq. ( p.la]) gives the retarded monopole moment, and ( p.lb| ) the retarded 2^-pole moment. 
Here, 

^i=yi-Zi{u), (2.2) 

and on P, 

e= = (e,k-e) , (2.3) 



or equivalently, 



kiC = . (2.4) 



The domain of integration P is the null hyperplane which contains the null geodesic from 
the field point x to C, and rfS is the invariant volume element on P, obeying 
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t^S^^y, (2.5) 

and it can be shown that 

dE^i = a-kkidT, . (2.6) 

Clearly, Tij,k^ j^^ is symmetric on all its indices after the colon, and by its construction, 
obeys 

fc'^%^..../„ = 0. (2.7) 

Integrals of the form (|2.1b|) arise naturally in this context in the following way. Points 

in the domain of integration of ( p. .181 ) lie on the past null cone with vertex x\ and so 

?/° = x° - |x - y| 

= + rk^ — |x — y| 

= ^° + k-e-a, (2.8) 

where 

|x-y|=r-k-^ + a, (2.9) 

and 

oo 

a = Y.hnr-^. (2.10) 

n=l 

The coefficients 6„ may be obtained from a power series expansion of |x — y| = |rk — ^|. 
Then if cr = 0, 

I/O = 2° + k ■ e 

= / + ^k.e, (2.11) 



so that according to (|2.4| ), lies on P. 
Introducing the notation used in GK, 



[F] = F(z° + ^k-e,y), 
it follows by applying the chain rule that 

[F\a = [FA + YA^A. (2.12) 

where [F] ^ means 'put = z° + and then differentiate w.r.t. y°'^ while [F^] means 

'differentiate w.r.t. and then put = + 

The conservation equation for the source ( |1.5| ) is used to obtain a sequence of conservation 
equations for the retarded moments as follows. In ( p.l2| ), we take F to be respectively 
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F = T^^ and F = T^^'^^i ■ ■ integrate over the intersection of the null plane P with the 
world tube W and use (|1.5| ) and the strong vanishing on the boundary to obtain 



kWTi, = , 



{ki:k2...k„) "I" ^"''^j:(fc2...fcn'^fcl) 

} = 0, 



(2.13a) 



(2.13b) 



To obtain the derivatives of the retarded moments, the procedure given by Hogan and 
Ellis (HE) [0 is followed. The calculation for the monopole moment is given in full detail; 
derivatives of the higher moments are obtained in a similar way. From (|2.1a|) and (|2.6| ), it 
follows that 



ij,k 



d 



Oj' kTijk}^ -\- 



d 



The latter term is given by 
d 



T- 



dx'^ 



T- 



-kk + -(^fc - (k • ^)kk) \ dsy . 



To write this in terms of the retarded moments, note first that 



In general. 



and using ( p^) , 



dF 

du 



Combining these results, one obtains 

^ ^ [T,,] = -kkD%, + ^{D%j.,u + v'kkDTija) 

and so 

%j,k = —%jkk + -{D%j:k + v^kkDTij,i) . 
Similarly, it is found that 

1 

%j:k,l - 

I 
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%j:kki + -{Vk{Tij.,i + v'^kiTij.m) + DTij-M + V^'kiDTij-km) , 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19a) 



(2.19b) 



and 



ij:kl,m 



+DTij.Mm + k^v^DTij-Mp) . (2.19c) 



Next, it is demonstrated how ( |1.18|) is evaluated. Again, the procedure of HE is followed. 
A power series expansion at infinity in the variable yields 



T.,(/,y) = E 



-a 



n=0 



T- 



where the subscript {n) means differentiate n times with respect to y^. Writing 

1 



-(J 



n\ |x — y| 



n = 0,l,2,.... 



( ILT^ ) becomes (in the IRF) 
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n=0 



T- 

J- 7 



«i>(n) 



dzy.. 



(2.20) 



(2.21) 



This can be written as a power series in r~^, with the coefficients identified as IRF values 
of certain combinations of the retarded moments. For illustrative purposes, the calculation 
for the first two coefficients is given. 

With the hn defined by ( p. 101) , it is found that 



whence 



1 1 



1 1,2 



k-^ + 0(r-3) 



H2 = 0(r-^) . 



(2.22) 

(2.23) 

(2.24) 
(2.25) 



Thus 
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I {4 1 [Td d,y 



+ - / 2k-e[T.,]-(|e|^-(k-0') T^UD dsy\ + 0{r 



The terms in the integrand here obey 

[Tij] dsy = [Tij] dZ = Tij , 

(ier-(k-o') = ier-(eY=e.f , 



kV 



(2.26) 

(2.27) 

(2.28) 
(2.29) 
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and so 



(ler-(k-e)^ 



ij:k ) 



(2.30) 



where ( p. 161 ), ( |2.17| ) and d[$,^] /du = have been used. Then combining these results, one 
obtains 
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ij:k 



+ 0{r- 



(2.31) 



Using the same procedure, the next term in the power series solution may be calculated 
to give 



r 



1 



'='^ + 0(r-^). (2.32) 

To end this section, it is pointed out that no assumptions, other than strong vanishing 
at the boundary, have been made about the source in order to obtain ( |2.32| ). In particular, 
the conservation equations (|2.13| ) have not been imposed, nor has a restriction to a source 
possessing only a finite number of non-zero retarded multipole moments been made. Finally, 
it should be pointed out that (|1.5|) is required to obtain the solution (|1.18|) , so that the 
comments below regarding withholding the conservation equations only strictly apply to 
quantities 7*^^ formally defined by retarded integrals of the form (|1.18|) . This will not affect 
the interpretation of the role played by the conservation equation (BTTBa). 



III. STRUCTURE OF THE CURVATURE TENSOR 

Having seen in the previous section how to evaluate the solution 7*^^ of (|1.2|), the curva- 
ture tensor associated with this solution is now evaluated and its structure examined. Sachs 



14| established the peeling-off theorem for linear fields with sources confined to a time-like 
world line; the same can now be done for the more realistic case of fields due to extended 
bound sources (see |^). This is interesting in its own right, but in addition will show up the 
important role played by the conservation equations in the determination of the peeling-off 
behaviour. 

Thus the main object of this section is to calculate the linearized curvature tensor (|1.6|) 



using the solution ( p. 32 ) correct to 0{r '^). Clearly, this would be an extremely lengthy 



procedure, so in order to establish the main results with the minimum fuss, the second 
derivatives of 7*^^ are calculated and then the NP Weyl tensor components are evaluated. 
This will give the following limited peeling-off result; 



r 



(0) 

^3 = ^ + 0(r-3) , (3.1b) 



(0) 

^2 = ^ + 0(r-^), (3.1c) 
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^i=0(r-^), (3.1d) 
^o = 0(r-^), (3.1e) 



where the = — 4 are calculated from the linearized curvature tensor ( |1.6| ). See 

e.g. [H for the definitions. These components will be evaluated on the null tetrad NT = 
{k\ n\ s\ s*} where 

1 1 dki _ 1 dki 

and so the only non- vanishing inner products are 

hn' = -1, Sis' = 1 . (3.3) 

The formulae of §1 and the expressions above for the first derivates of the retarded 
moments are used to obtain the following expressions. (Note the use of (|1.13|) in these 
calculations, and p = —a-kp, which follows immediately from ( |1.11|) .) 



tij^k = -%jh + J(4 + v'kk)D%.,i , (3.4) 
%j:k,i = -%r.kh + \5]^ + v"'ki){vkDTij.,m + OfcT^j:™ + D^Tij-km) , (3.5) 

Tij-.m, I — Tij:m ki ~\~ ~\~ V^ki)(2v DTij ^np ~\~ 2ci Tij-rnp ~\~ D Tij-m p) , (3-6) 

= -f^r.m"'kl + ^{Sf + v^h){2D\v'^%,.,^p) + D^%,.,rrp) ■ (3.7) 

Using these and the expressions for the first derivatives given in §2, the second derivative 
of the (trace reversed) metric perturbation can be calculated. The resulting lengthy expres- 
sions are given in Appendix A. Using these, the complete tensorial expression for Lijki may 
be given, correct to 0(r~^). However, the main results may be expressed more succinctly in 
the NP formalism. 

Foremost of these results is the verification of ( p.l|) , which is done using the expressions 
(|A6|)-( |AT0D of Appendix A, ( |1.6| ) and the relation ( p..4| ) between and Yij- the general 
case, it is found that 

*f = -2{D + a-k)DTijs'P , (3.8) 



= 4a ■ sk'n^Tij + a ■ sT; - 2s\a^ + a-kv^)Ti 



-2s\D + a-k)DTj^."' + s'{D + a-k)DT\.,j . (3.9) 

An examination of Appendix A will convince the reader that ^^2°"* is too long to merit 
inclusion. For simple illustrative purposes it may be calculated for the case of a monopole 
particle, that is for a source for which 

Tij ^ 0, %j-M...kr, = n > 1. 

In this case, it is found that 
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777 

v^, = -^, (3.10) 

where m is a constant, which is the usual resuh. 

It should be noted that in deriving eqs.( |3.1| ) and the coefficients above, the conservation 
equations (|2.13|) have been used repeatedly. A more careful examination of the leading order 
term of the field demonstrates what role these equations play. It is found that 



■'ijkl 



N,,kir-' + 0{r-^) , (3.11) 



where 



and 



Nijki = Niikjkk + Njkkiki - Nikkjk - Njikikk , (3-12) 



= {D + a-k)DTij - ^7]^j{D + a-k)D%n"' . (3.13) 
Imposing the conservation equation ( p.l3a| ), it is found that 

Nij^ = -i(D + a-k)D%rki , (3.14) 

and hence 

Ni,Mk' = 0. (3.15) 

Thus, to leading order the field is Petrov type N, with degenerate principal null direction /c*. 
This vector has geodesic, twist-free, shear-free and expanding integral curves, the null rays 
of the radiation field. This will allow in §4 a discussion of the geometric optics approximation 
to the field, by considering the field 

^ijkl = — ^ , (3.16) 

and working along the same lines as §4 of 0]. 

It should be pointed out that in general, i.e. without imposing the conservation equations 
( ^.13| ), the leading term in the field ( |3 . 1 1| ) is not type N. Indeed 

Ni.kik' = D{k'D%i)k,kk - D{k'D%i)hkk, 

which will be non-zero in general. However, it is straightforward to show that 

Niji^kkm]k' = Q , (3.17) 

so that the leading order term is Petrov type III, with degenerate principal null direction fc*, 
which type of field is still characteristic of gravitational radiation. Thus it is seen how a fun- 
damental property of the source, the monopole conservation equation, has consequences for 
the dynamics of the field. These consequences of conservation will also manifest themselves 
when the flux of energy-momentum of the field is considered. 
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IV. THE GEOMETRIC OPTICS APPROXIMATION 



In this section, the geometric optics approximation to the Lorentz covariant gravitational 
field is discussed. Two related aspects of this problem are considered. Firstly, it is shown 
that there exists a linearized curvature tensor which may be considered as an approximate 
solution to the field equations, and which has the geometric structure of a radiation field. 
Secondly, the energy momentum of the full field is examined, using the Landau-Lifschitz 
(LL) pseudotensor This is shown to split into a geometric optis parts, which gets radiated 
away to JT""*" along null geodesies, and a remainder term which in the way set out below, 
remains bound to the source. 

The LL pseudotensor appropriate to linearized gravity is M 



±ik ^ ^ r)-pnik-pp -pni ppfc pm -ppk 



which is conserved, 



/cy-pnrn pp pnm, pp pnm pp \ ik 

m-L np ^ pJ- mn ^ mp) '/ 

I pfei ppn I pfcn ppi pfcm ppi pfei ppm 

"1"-'- p^ n ~T~ ^ p ^ p^ m ^ p 

ipifc ppn I pin ppfc pirn ppfc pifc ppm 

"rJ- pJ- n ~r J- n-L p J- pJ- m J- p 

_i_pim pfc " p«m pfcn \ (AT^ 



f\, = 0. (4.2) 



Obviously, this is not constructed from the field Lijki, but from the 'field strengths', T'^jk- So 
unlike the situation in electromagnetic theory, where the geometric optics part of the energy- 
momentum tensor of the field is the energy-momentum tensor of the geometric optics part 
of the field, the geometric optics approximations must be defined separately. The relation 
between the two will be clear. A different approach to the geometric optics approxiation in 
linearized gravity involves a WKB analysis of the field equation ( |1.2| ) 

The first part of the discussion involves an analysis of the leading order term in the field, 

^ijkl = ~^ijkl , (4.3) 



where Nijki is defined in ( |3.12D and ( p.l3| ) 



To show that ([4.3|) may be considered an approximate solution of the linearized field 
equations, for large values of r, it must be shown to be (approximately) a linearized Rie- 
mann tensor with vanishing Ricci tensor. Away from the source, the linearized Ricci tensor 
constructed from (|4.3| ) is 

n,,=C,,k' =0. (4.4) 

Now Cijki clearly has the symmetries of a Riemann tensor, so it remains to show that it 
obeys (approximately) the linearized Bianchi identities, 

Lij[ki,ni] = . (4.5) 
Using some results given in §3, one can calculate 
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ij,k 



-2{D - a-k){D + a-k)DTijkk 
2 

+- {(4 + v^h){D + 2a-k){D + a-k)DTij.,i + 3{ak + a-k{vk - kk))fij 



+ {bk + b ■ vkk + Qa-kak + {b-k + 6a-k'^){vk - kk))Tij^ . 

This is the central equation used to obtain 

1 



ij[kl,m] 



^Cijl^klkm] ; 



where 



Cijki — 8B[ij][ki] , 

Bijki = {{D + 2a-k){D + a-k)D%l.i + 8(0^ + a-kvi)%l 

+ {bi + 6a-kai + {b-k + 6a-k'^)vi)T*^} kj + (D + a-k)DT*{riji + 2k ^vi) 

A rest frame calculation along the lines of §2 shows that 



ATI* 



(4.6) 
(4.7) 

(4.8) 
(4.9) 

(4.10) 



Bijki = kj [t;;,(3)0] 4z/ + {riii - '2kj5f) 
The domain of integration here is the sphere 

ICI < Ku), 

where h is the 'radius' of the source given in the introduction. Also, since Tij{x) is, in each 
integral, compactly supported, there exist scalars Mn{u), n = 0, 1, 2, . . . such that 



|7;;.(„)| <M„(m), n = 0,l,2,... 



Using this fact, the following bound can be obtained for ( [4.10|) ; 



\B,,ki\ < 47r/i3M2 + ^tt/i^Ms 



(4.11) 



(4.12) 



which leads to a similar bound for Cijki- 

Thus equations ( |4.4| ), ( |4.7| ) and ( |4.12|) show that Cijki is an approximate solution of the 
vacuum linearized Einstein equations. The accuracy of the approximation is determined by 
the bounds M2 and M3. 

It has already been shown that Cijki is Petrov type N, with fc* as repeated p.n.d. There- 
fore, the field picks out a unique null direction whose integral curves form a shear- free, 
twist-free and expanding congruence of null geodesies. This ray geometry is characteristic 
of a geometric optics field, and so Cijki is referred to as the geometric optics approximation 
to the linearized gravitational field. One usually thinks of the field as 'propagating in the 
direction fc", and so it is natural now to consider the question of energy- momentum transfer. 



For this, the LL pseudotensor (4J^) is used. As pointed out above, this is not a functional 
of the field (curvature tensor), and so a seperate but related geometric optics approximation 
is assumed. A detailed examination of the first two terms (coefficients of and r~^) of t^^ 
is required, to obtain which one can calculate 
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4 4 
l*ij,k = — DTijkk + — {Tij{vk - kk) + DTij,k 



+ {v\2D + a-k)%j.,i + a%j.,i + ^{D + a-k)D%j./)kk^ + 0{r-^) , (4.13) 



which leads to 



r,,., = -l{DT*kj + DT*k, - Dr;,ki) 

+0{r-^) , (4.14) 

where 

a,j = ^(D + a-k){D + v\W + a-k)^., + a'T;*, - r*. . (4.15) 
Trace reversals are taken over indices before the colon. One then finds 

-{\iD + a-k)DTZl + v\D + \a-k)TZi + ^a'T^, - ^T^)^, . (4.16) 

Using this, it is found that 

r'^^ = 7**^\, = , (4.17) 
r^, = '^-DT,% - ^{DT;^. + T/vj - a\kj) , (4.18) 

^^Jkk'' = + 7;VJ^ + b[n]) + Oir-"") , (4.19) 

2 2 f 1 1 

V^jk^' = -DT^kjkk - iTjk + -VjkT/ + -T/{vjkk + Vkkj) 

\DTZ,kk + \DrZkk, + {DTki! + Tkiv^)k, + {DT,i! + T,iv^)kk 

+{\{D + a-k)DTZl + v\2D + a-fc)r- + a'T- - TZ)k,kk] + 0(r-=^) . (4.20) 

where 

h^i = (r>, + DTZi - 2DTimr - '2Ti^v^)kj , (4.21a) 

which gives 

bjik^ = bjik' = T^kj . (4.21b) 
It will be convenient to write 



oo j-ij 

t'' - E Sir ' (4-22) 



T 

n=l ' 
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where 



and to define 



where 



dr 



n > 1, 



rad^ ~l~ rem^ 



dU 



{l)Uj 



(4.23) 



(4.24) 



(4.25) 



The 'rad' and 'rem' here stand for 'radiation' and 'remainder' respectively. Using ([4. 1|) and 
( CT ) - (iHI), it is found that 



rad 



27rr2 



(4.26) 



which is the usual form for the energy-momentum of a geometric optics or radiation field. 
It is straightforward to show that this obeys 



rad^ . 7 



and since the total pseudotensor is conserved, the remainder is also conserved; 

0. 



rem j 



OO 



\n=2 



+1 



(4.27) 



(4.28) 



Now a long calculation using ( [4. 1|) and ( |4.14|) - ( [4.21|) yields the following important result, 

f^kj = 0{r-^) , (4.29) 

which from ( [4. 261 ) is equivalent to 

(2)f ^fc, = . (4.30) 

Hogan and Ellis showed that the energy-momentum tensor of the electromagnetic field due 
to a bound source obeys a set of equations similar to ( [4.26 ) - ( 4.30 ) above. Thus a theorem 
of theirs can be applied to the linearized gravitational field: 



Theorem 1 Let K'^" = U'^k^ — W^k", where 



U 



ik 



(n) 



ik 



1 °° 1 

^^^(2) ^" \n — 1 {n — l){n — 2) 



Then one can write 



,k 



n=2 ' 



fij 

1 ^ 



(4.31) 
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Proof. See g, p.202. 

To examine the transfer of energy-momentum by the field, one considers the flux of 
4-momentum across the following fundamental surfaces. This follows exactly the procedure 
of HE for the electromagnetic case. 

Consider the 4- volume V of Minkowskian space time bounded by two time-like 3-surfaces 
r = ri and r = r2 > ri and by the null cones u = Ui and u = U2 > ui. ri is taken to be 
large enough so that we are outside the source, and so away from the world-line C given by 
r = 0. Thus f^^ is non-singular in V . The fluxes of 4-momentum across the boundaries of V 
(in the directions of increasing r, u on the appropriate surfaces) are given by 



{A)pi _ ^2 / ^ij^ , ^ ^4 32) 

J ui J 

across r = rA_, {A = 1,2) and 

(^)g* = - J^' dr J f^kj du , (4.33) 

across u = ua- In these integrals, du is the area element on the unit 2-sphere. The 
conservation equation = gives the conservation law 

Wpi +(1) Qi = (2)pi _^(2) Qi ^ ^4^34) 

which says that the total 4-momentum entering V is equal to the total 4-momentum leaving 
V. 

Using the decomposition of the energy-momentum pseudotensor ([4.24|) , one can write 

P' = r.dP' + rernP\ (4.35) 

on r =constant, and 

Q'=radg^+rcmQS (4.36) 

on u =constant, where the first terms on the right hand sides are obtained by using rad^*"* 
in the definitions ( [4.32| ) and ( [4.33| ), and the second terms on the right hand side by using 

fij 
rem'' • 

From ( [4. 261 ), one finds that 

1 /•«2 /• . 1 



radP^ = ^ £ J {DT^nDT"^" - -(I)Tr)2)A;^ duj , (4.37) 



and 



radQ* = . (4.38) 

Noticing that rad-P* is independent of r, these results give trivially the conservation law 

lUr +l'lQ' = i^r +^^,Q\ (4.39) 

which also follows from ( [4.27| ). This also means that rad-P* may be used to calculate the 
flux of 4-momentum at J'~^ (r oo, u finite). Indeed, since fc* is future pointing and the 
coefficient 
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m\2 



> 



(see below), one sees that rad-P* is future pointing, and so this vector is referred to as the 
total radiated ^-i^omentum of the field across r = constant in the retarded time U2 — ui. 

At this stage, the question of outgoing radiation conditions may be addressed. Traut- 
man has expressed Sommerfeld's outgoing radiation conditions as follows ||20|,|3[| : there exist 
coordinate systems and a tensor fij = 0{r'^) such that 



9ij = Vij + 0{r ^) , 



9ij,k 



f:^k^ = 0{r-^) , 

with /j* obtained from (pT^). 

The energy-momentum pseudotensor calculated from ( |4.4CI|) turns out to obey 



(4.40a) 
(4.40b) 
(4.40c) 



f3 



327r 



J mnJ 



if 



*m\ 2 



k'y + 0(r- 



(4.41) 



and the coefficient of here is positive. Indeed, given any symmetric second order tensor 
Aij on Ai obeying Aijk^ = 0, one can write 

= cakikj -\- P(^kiSj -\- kjSi) -\- (3{kiSj -\- kjSi) 
+7SiSj + ^SiSj + 5{siSj + SjSi) , 

where a and 5 are real valued and the null tetrad of §3 has been used. Then A\ = 25 and 

AijA'^ = 277 so that 



Aij A 



n2 



277 



which is positive. This is applied to f*j, using ( [4.40c| ), to show that coefficient of /c* in 
the integrand in (f4.37|) is positive, so that the integral leads to a positive outward flux of 
4-momentum of the fleld. 

Now from ( |4.13| ), it is seen that ( [4.40a| ) and ( [1.40b| ) will be satisfled for the fleld by taking 



fij = -{4D%,-27]i,DTk'')r-K 



(4.42) 



Then the key equation ( ^.40c| ) which is needed for the procedure above is 

kWTij = , 



(4.43) 



which is exactly ( |2.13a| ), the monopole conservation equation of §2. 

This may be compared with the result obtained by Hogan for a moving monopole particle 
[HI. Consider then a source for which 



From (|2.13b|) with 



n 



1, one obtains 



ij'.k\ . . .kfi 



0,n > 1 



(4.44) 
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%j = mViVj , 



for some scalar m. Then the outgoing radiation conditions, i.e. 

(lijk^^ = —mvi — mai = , 



give 



(4.45) 



(4.46) 



from which (contract with f*, then a*, using f jO* = 0) it is found rh = and a* = 0. This is 
exactly the result obtained by Hogan; the outgoing radiation conditions force the acceleration 
of the particle to be zero. What is seen here is how this relates to the conservation of energy 
momentum of the source, i.e. of the particle. In general, this shows how a fundamental 
property of the source, its conservation, has important consequences for the dynamics of the 
field. 

Next consider the remainder terms in (f4.35|) and ([4.36|) . Using ([4.31|) and Stokes' theo- 
rem, one can write 



K^^^kjVk duo 



U=U2 



U=Ui 



and 



+ 



.2 J X'^kf^^^^ 



'• f K'^^kjVk du 

Using the definitions of the theorem and (|4.29 ), one can write 



K^^%v,du; = j:-^ 

n=3 ' 



fn-2) 



where 



/;(u) = / ^4'^kjduj, n = 3,4,5, 



This allows one to write 



/n(«2) - fiiui) 



n=3 ^1 



n-2i 



(2) pi = V 
rerir / ^ 



n-2) 



n=3 



and 



71=3 

oo 



(2) ni 
rem 



n=3 



[n-2) 



1 



"2 / 
1 \ 



n-2 



^"-2 I ' 
'2 / 



which lead to the expected conservation law, 

(1) pi , (1) Qi ^ (2) pi , (2) Qi 
renr ^rem^ rerir ^reiiiT' 



(4.47) 



(4.48) 



(4.49) 

(4.50) 

(4.51a) 
(4.51b) 



(4.52a) 
(4.52b) 

(4.53) 
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Unlike the case for rad-f", the fluxes in ([4.51|) depend on r. Eq. ( [4.52| ) indicates a leakage of 



4-momentum across null cones, so that the remainder 4-momentum is not purely radiative. 

One can make the same observations which Hogan and Ellis made for the electromagnetic 
field. The total flux of 4-momentum across r = constant in a proper time U2 — Ui is 



rem-* 



= pKu2) - pKu,) , (4.54) 
with 



00 



This vanishes in the limit r — 00 and so is non-radiative in character (no flux at J^~^). This 
qualifies the characterization of rad-P* as being the total radiated 4-momentum. Eg. ( [4.54 ) 



may be said to describe 'particle-like behaviour', in that the flux of 4-momentum between 
u = Ui and m = M2 is simply the difference in the values of the '4-momentum ' p^{u) at the 
two times. For this reason we refer to p*(-u) as the bound 4-momentum of the source 0,|18 



Thus Teitelboim's idea of the bound 4-momentum of the electromagnetic field [0 may be 
generalised to (at least) the linearized gravitational field. 

The splitting of the flux of 4-momentum across r = constant surfaces defines the ge- 
ometric optics approximation for the energy-momentum of the field. The relation to the 
geometric optics part of the field (curvature tensor) itself is clear; only terms which are 
involved in the construction of Cijki are involved in the construction of rad'^*-'- Thus asso- 
ciated with Cijki is the total radiated 4-momentum, which detaches itself from the source 
and radiates away to JT""*" at the speed of light. The energy-momentum which is left behind 
remains bound to the source. 



V. ARBITRARY DEPENDENCE ON A NULL COORDINATE 

In this section, results are obtained which can be used to determine the relation between 
the retarded multipole moments of the source and the multipole moments of the field. This 
will involve determining a decomposition of the field into arbitrary functions of u (which is 
a null coordinate), and known functions of angle. As in the electromagnetic case, radiative 
linearized gravitational fields must contain arbitrary functions of a null coordinate in order to 
convey information ||2l[]. The decomposition obtained will allow a rewriting of the solution 



7jj in a form which allows direct comparison with other approaches to linearized gravity 



In order to keep the calculations to a reasonable length, attention is restricted to a 
source for which only the retarded pole, dipole and quadrupole moments are non-zero. 
Such a source will suffice to display generic behaviour. Thus 

'TirM...k„ = , n > 3 . (5.1) 

The decomposition of the retarded pole, dipole and quadrupole moments is obtained by 
analysing the conservation equations ( |2.13| ). 
Defining the tensors 



19 



ITT'i — Tijk'^ 5 'm'i:ki...k„ — ^"'^j:fci...fc„ ; (5.2 

the first few conservation laws may be rewritten as 

rhj = , (5.3 

mi.,j + Tij + rriiVj = , (5.4 

mi-.jk + + Tik:j + rrii.jVk + mi-kVj = , (5.5 

rhi-jki + %j:ki + %k:ji + %i:jk + mi.jkVi + nii-MVj + mi-jjVk = . (5.6 
Tlie first of tliese may be rewritten as 

1 1 

Tij = --{rhi.,j + rn.j.,i) - -{rriiVj + nijVi) . (5.7 

From (|5.4|), one can deduce tliat tliere exists a tensor Qij^ such that 

Tij;k ~ Tf^i-.j'Vk '^T^i-.jk Qijk i (^-^ 

where Qijk obeys 

Qijk Qikj 5 (5.9 

Qijkk = Qikjk — ■ (5.10 

Similarly, using (|5.5| ) along with (|0|), one finds that there exists a tensor Qijki such that 



%j:ki = --{vkmj;ii + virrij-ik) + Qjiki , (5.11 



and Qijki obeys 



Qijkl — Qijlk 5 Qijki + Qiklj + Qiljk — , (5.12 

Qijklk = — — nT'i:jk^ Qijklk'^ = ^i:kl ^ (5.13 

Qijkik^k^ = Qijki^k^ = Qijklk'' k^ = . (5.14 

The procedure runs as follows: using the formulas for the derivatives of the retarded 
moments, (^-derivatives of Qijk and Qijki are obtained. These give a pair of first order 
simultaneous differential equations, which are in an integrable form. The 'constants of 
integration' yield the required arbitrary functions of u. 

To begin, recall that 
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d dk^ d 
r- 



Then using the derivatives (|2.19|) one finds 

dmi:jk _ ,^ dk^ 



and so 



dC 



dC 



dk^ 



Using (|2.19b|) and ( p.l9| ) one obtains respectively 



(5.15) 



(5.16) 



and 



dm 



dk^ 



Using these derivatives and the definition ( ^.8|) of Qijk gives 



(5.17) 



(5.18) 



dQ 



ijk 



ld_ 

2du 



dk^ 



(5,19) 



Similarly, one can use ( p.l9c|) and ( |2.19| ) in the definition ( ^.111) (remember = v^{u)) to 
obtain 



dQj 



jikl 



1 dk"^ 

2 {i^jidni Tjl:irri)'^k ~l~ i^ji.km '^jk:im)^l\ 



(5.20) 



Now from ( ^.19] ) one can deduce the existence of tensors Rijk and Qijk{u), both having the 
same symmetries as Qijk, and Qijk depending only on u, such that 



Qijk Riik Q 



and 



dRijk 



ijk 



hr -r )^ 

^{■^ij-.km ^ik:jm) 



Comparing this with (|5.20| ), one can write 

Qjikl = —RjilVk — RjikVl + Qjikl{u) , 

where Qjiki has the same symmetries as Qjiki and depends only on u. 
Eq. (|5.22| ) is solved as follows. Using ( ^.23| ) and the first of ( ^.13|) , 



(5.21) 



(5.22) 



(5.23) 
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which on using ( |5.23| ) and ( |5.11| ) gives 

Tij-.kl = '^Rijmh VkVl + Qijkmk Vl + Qijlmk Vk + Qijkl ■ 

Then 

{T\j:kl T\k:jl) {Qijkl H" Qijlmk Vk Qikjl Qiklmk Vj) 

d 



(5.24) 



(5.25) 



^ (Q^Jlm{Sk + ^k^Vk)k' - Q^kiU^T + ^k^v,)k' 



Comparing this with ( 5.22| ), one can integrate and give 

R-ijk 



Q^M^k + lk^Vk)k' - Q^kU^f + h^V,)k' \ . 



(5.26) 



(5.27) 



The constant of integration which arises here can be absorbed into Qijk and Qijki, and can 
thus be set equal to zero without changing Qij^. or Qijki- From ( |5.21| ) and (|5.23| ), 



Q 



ijk 



QvU^T + lk^Vk)k' - Q^kim{5f + \k^v,)k'^ + 



and 



'2du 



Qj^kl = -Vl Q,ipm{5^ + -k'^'VuW' - Q.kpmK + -k"'Vl)kP 



(5.28) 



+Qjiki{u) ■ 



(5.29) 



The tensors Qijk and Qijki will be specified by the source via these last two equations and 
( ^.8|) and ( p.ll| ). The linear transformation Q = Q{Q) implicit in ( ^.29| ) is not invertible, 
that is Qijkl = does not imply Qijki = 0, but one will certainly be able to determine the 
non-pure gauge parts of Qijki from this equation in terms of Qijki, and so in terms of the 
source. This will be seen more explicitly when the solution 'jij is written in terms of these 
tensors. 

To do this, the combinations of retarded moments which arise in the solution ( 2.32| ) are 
obtained in terms of rrii, Qijk and Qijki- For a pole-dipole-quadrupole source described by 
( ^yri) , this solution can be written exactly as 



7 



^7-, - \{2%,..kv'' + D7-,:/) + -(2%,.,kiv%' - %,.,k') . 



From (|5.29| ) and (|5.11|) , one has 

Tij-.kl — Qijpmk^k VkVl + Qijlpk^Vk + Qijkpk^Vl + Qijkl ; 

so that 



(5.30) 



(5.31) 
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%j:k' — —Qijkik^k^ + '^Qijkiv^k^ + Q\jk -I (5.32) 

and 

%jkiv^v^ = Qijkik^k^ - 2QijipV^v'P + Qijkiv^v^ . (5.33) 
Using (15^281 ) and (15:211) , one can write 



%j.,k = -2Rijik^Vk - a-kRijik^Vk - Rijik^k + k^DQijki - Qijik^v^ - Qijk , (5.34) 
and so 

Tij.kv'' = 2Rijik^ + a-kRijik^ + {Qijkik^yv^ + Qijik^ - Qij^v^ , (5.35) 



and from (|5.10|) 



rrii-.k = Rijkk'' + Qijkk^ ■ (5.36) 
Then a straightforward calculation gives 



= iQijkik'^k' + ^a-kQijkik'^k' + ^{b-k + 3{a-kf)Qijkik''k' 

Qijik'' — a-kQijik^ — rriiVj , (5.37a) 



+Qijki{a-k{7]''^ + Qv^k^ - Qk^k^) + 2a}'k^) 

+2Qiji{k' - v') , (5.37b) 



2%,.Mv\' - %r.k^ = Qi.kii-v"' + ^k'k' - Qv'k' + 2v''v') . (5.37c) 



Thus from these expressions and ( [5.30| ), it is seen that one can express the solution in 
terms of arbitrary functions of u, namely rrii (which is in fact constant), Qijk{u) and Qijki{u). 
One can go further than this, and show that these tensors act as potentials for the solution. 
A direct calculation shows that 

* — (Q ki + a-kQ'^'ki) + -Q'^\vi-ki), (5.38) 



r 



{Q +3a-kQ +{b-k + 3{a-kf)Q'^^^)kkki 



1 / \ ijkl 

+— Q {-rjki - Avkk + 3kkki) 



+ Q'^''\-a-k{7]ki + 6vkki - 6kkki) - 2akki)^ 

+ 4 iQ'''''i-Vki - Gvkki + Skkk + 2vkVi)) . (5.39) 
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Comparing these with ( |5.37|) and ( ^.301) above, it has been shown that 



^.«_,^^jQ![Y,+ f^),.. (5.40) 

r \ ^ J V / 

This form of the solution can be used to show that Q*"'^ ^ind Q*"''^' contain the same infor- 
mation as Q^^^ and Q'^^^K To see this, notice that 

/ r>iik\ 

(5.41) 
and 






. ^ / V ^ / 



r 



M 




(5.42) 



where {v'^ lr)^i = and Rijk = Rijkiu, (, () have been used. Thus 



2(«!!),W«!^),,.2(«^),+ f«!^)„. (5.43) 

\ r J \ ^ J ' \ ^ / \ ^ / 

and so ( p.4(J| ) may equally be written without carats on the tensors Q''^^ and Q''^^K This 
validates the comments after Eg. ( ^.291) ; the kernel of the transformation only yields terms 



which do not contribute to 7**-^, i.e. pure gauge terms. 

Using the tensors Q^^^ and Q^^^'', one could calculate what Janis and Newman |Q call the 
multipole moments of the field (the coefficients of inverse powers of r of the NP components 
of the curvature tensor). These would be given in terms of f*, the Q tensors and known 
functions of angle (from the null tetrad of §2). This lengthy calculation would result in 
combinations of our arbitrary functions of u and such spin-weighted spherical harmonics as 
are predicted by the general theory [p!0| , pJ]] . 



The solution ( |5.4(j| ) is equivalent to those given by Sachs and Bergmann [T^l and Pirani 
||T3| for 'multipole particles', but has the advantage of being derived from a realistic ex- 
tended source rather than one confined to a time-like world-line. These authors make use 
of symmetric trace-free tensors which are functions of the retarded time, and their solutions 
contain infinite series of such tensors, corresponding to a particle possessing an arbitrary 
number of multipole moments. It is clear then that for a general source for which arbitrar- 
ily many of the retarded multipole moments ( p.l|) are non-zero, the solution (|5.40|) would 
contain an infinite series of tensors Q'^i'"''"(u) which are arbitrary functions of u. Each of 
these would be present at the level in jij, and would lead to an infinite series of spin- 
weighted spherical harmonics appearing in the fully decomposed expression for T^. As seen 
in the previous section, in some situations, and especially those involving a discussion of the 
radiation field, it is advantageous to use the retarded moments (|2.1| ). 
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VI. CONCLUSIONS 



The similarity between the electromagnetic and linearized gravitational fields has been 
exploited in this paper to derive some results about the latter using the Goldberg-Kerr 
approach to bound source fields. The results refiect those obtained by Hogan and Ellis 0] 
for the electromagnetic field. Well established results (e.g. peeling-off) have been confirmed, 
relying on the use of the retarded multipole moments of the source. The advantage of 
using these can be seen in, for example. Eg. ( [4. 371 ), where the total radiated 4-momentum 
is determined by a single integral over the source, the retarded monopole moment. The 
results of §5 show how this relates to analogous formulae of |TP[, which involve infinite 
series of moments of the source. In addition, use of the retarded moments of the source to 
describe the exterior field has made the consequences of the conservation of the source easy 
to identify. 

It comes as no surprise that a bound source linearized gravitational field possesses a 
bound 4-momentum. The obvious question is: does this notion extend to the full non-linear 
gravitational field in the case of asymptotic flatness? The Landau-Lifschitz pseudotensor can 
be used to construct the 4-momentum for such fields (and the calculations would be carried 
out in an asymptotically Minkowskian coordinate system, c.f. §20.2 of 0), and indeed this 
4-momentum conincides with the covariant (Bondi-Sachs) constructions [Q. Hence one 
might expect that some or all of such fields do possess a bound 4-momentum, which could 
be constructed in a manner similar to that above. 
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is determined by a single integral over the source, the retarded monopole moment. The 
results of §5 show how this relates to analogous formulae of |jl9|, which involve infinite 



series of moments of the source. In addition, use of the retarded moments of the source to 
describe the exterior field has made the consequences of the conservation of the source easy 
to identify. 

It comes as no surprise that a bound source linearized gravitational field possesses a 
bound 4-momentum. The obvious question is: does this notion extend to the full non-linear 
gravitational field in the case of asymptotic fiatness? The Landau-Lifschitz pseudotensor can 
be used to construct the 4-momentum for such fields (and the calculations would be carried 
out in an asymptotically Minkowskian coordinate system, c.f. §20.2 of 0), and indeed this 



4-momentum conincides with the covariant (Bondi-Sachs) constructions [0. Hence one 
might expect that some or all of such fields do possess a bound 4-momentum, which could 
be constructed in a manner similar to that above. 
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APPENDIX A: SECOND DERIVATIVES OF THE METRIC 



In order to obtain the second derivatives of the metric perturbation 7,^- correct to 0{r~ 



the solution (|2.32|) is written as 

■y', = ^ + ^ + ^ + 0(r-'), (Al) 

where the coefficients may be read off from ( 2.32|) . For convenience, the decomposition 

Bij = ttij + Pij + fiij , (A2) 

will be used, where 



a 



fJ'ij 



-4T- v"^ 

^■'■ly.m^ 1 
ij:m 1 



(A3) 
(A4) 
(A5) 



Then all the information required to construct the field (to 0(r~^)) is contained in the 
second derivatives of the terms Ajj, Bij, aij, (3ij and /ijj. These derivatives are obtained 
using results given in §§1-3. It is found that 



( ^ 

V r 



{{D + a-k)D%j)kkh 



--{2{D + a-k)D%,.,U5'(l + v^k(k)h) + D%,.,m{(^^ + a-kv^)hku 



DTij{r]ki + Avi^kki-) - ?>kkki) + a-kTij{2vi^kki) - ^kkh) + 2Tij{ai^kki))} 



+^|'7ii(-^?M + '^.VkVi - 6v(^kki) + Skkk) + v'^D%j.,ra{r]ki + '2,V(kki) - kkh) 
+2%,.,ma(^k{Sl^ + v'^h)) + 2{2D + a-k)%,U5^ + v^ki^i){vu) - k^)) 
+{D + a-k)D%^.,^,{5'^ + v^kkWi + ^'h)} , (A6) 



= -^{{D + 2a-k){D + a-k){v"'%,.,,^)} k^h 

4 r 

+2v"^%j.m{akki + aikk + a-k{r]ki + '^Vkh + Avikk - Skkk)) 
+ {~{D + 2a-k)Tij.,p + v'^iD + 2a-k)Dij.,^p + a'^Tij^mp) x 
{{6f + vPki)kk + {6l + vPkk)ki) 

+ {aP + a-kvP)i-%j..p + D%j.,^p)kkki} + ©(r"^) . (A7) 



,ki = -^{D + 2a-k){D + a-k) fij-m'^ kkh 

2 ■■ m 

+ — |Ty:m™ {riki + ^Vkh + ?,vikk - hkkki) 
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+2 tij-m"" {akk + aikk + a-k{r]ki + ^Vkh + Avikk - Skkk)) 
+{2(D + 2a-k)D(v'^%j..^p) + (D + 2a-A;)L>%:p^"^) x 
((5f + vPk)kk + (5l + vPk,)ki) 

+{aP + a-kvP){2v"'D%j..^j, + 2a^%^,^p + D'' %^,^^)kkki} + 0(r-^) . (A8) 



= - — {D + 2a-k){D + a-k){a-k%j.,rn'^)kkki 

+^{2(6-A; + 3a-e){2v^Tij.,mp + D%^.,pm'^){5l + v^k(i)kk) 

-\-2a- k(^2v"^DTij-jjip + 2a"^Tij.mp + D'^Tij.p^ ^ '^^^k{i)kk) 
+a-k{2v"'%,.,mp + DTij^prnm^^" + a-kvP)kkki} + 0(r-^) . (A9) 

= ^{D + 3a-k)(D + 2a-k)Cijkkki + ©(r"^) . (AlO) 
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